      SUBROUTINE INLOD(X1,Y1,X2,Y2,DG,NEAP,ID,J,alfa)

C
C ESTA SUBROTINA CALCULA O VALOR DOS ELEMENTOS DA DIAGONAL DA MATRIZ dG
C

      IMPLICIT DOUBLE PRECISION (A-H,O-Z)

      DOUBLE PRECISION L

      DIMENSION L(2)

      L(1)= (X2 - X1)
      L(2)= (Y2 - Y1)

      XL= DSQRT( L(1)**2 + L(2)**2)

      SIGN= 1.
      IF (J.EQ.NEAP) SIGN= -1.

      DG= -SIGN*(alfa+ DLOG( 2/XL ))*L(ID)/XL

      RETURN
      END